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Importance of high p t and nuclear physics for simulating UHECR 

Air Showers 

Jeff Allen*, Glennys R. Farrar* 

* Center for Cosmology and Particle Physics, New York University, USA 



Abstract. Observational evidence from Auger and 
earlier experiments shows a deficit of signal in a 
surface detector compared to predictions, which in- 
creases as a function of zenith angle, when the energy 
of the event is fixed by fluorescence measurements. 
We explore three potential explanations for this: the 
"Cronin effect" (growth of high-transverse momen- 
tum cross sections with nuclear size), the need for 
more particles at high transverse momentum in p—p 
collisions than currently predicted by high energy 
hadronic models used for air shower simulations, 
and the possibility that secondary interactions in the 
target air nucleus produces additional relatively soft 
pions not included in simulations. We report here 
on the differences between Pythia and QGSJet II, 
especially for high p t particles. The possible impact 
of these effects on the predicted surface array signal 
are also reported. 
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Fig. 1: The "Cronin effect" is the growth of a, defined 
in the text, with p± such that it significantly exceeds 
a = 2/3 CQ . 



I. Introduction 

A discrepancy between observations and simulations 
of ultra-high energy air showers has been recognized 
for many years, initially evident as an energy-scale 
discrepancy between air-fluorescence and surface array 
experiments. A discrepancy is also seen between the 
predicted and observed "attenuation curve", plotting 
some measure of the surface array signal as a function of 
angle, for a constant intensity cut (e.g., for the N highest 
energy events at the given angle). The discrepancy can 
be parameterized by a "muon rescaling" factor and 
an energy rescaling factor, En esc and the results of this 
are given using a number of different techniques in the 
presentation of A. Castellina for the Auger collaboration 
at ICRC 2009. So far, efforts to remove this discrepancy 
by modifying the p-p event generators (QGSJet, EPOS, 
Sybill,...) have not succeeded to fit both the Auger 
observations and other data. 

Here, we investigate three different examples of "ordi- 
nary physics" that are not included in present event gen- 
erators, which can increase the signal in a surface array: 
1) the "Cronin effect" 0], 2) a systematic deficiency 
in the production of high p± particles in "minimum- 
bias" event generators like QGSJet compared to event 
generators like Pythia which have been tuned for high 
P± El, 0, and 3) the possible production of an excess 
of low-energy pions in the target-fragmentation region 
due to secondary scattering in the target nucleus. 



II. Cronin Effect 

In the 1970's, it was discovered that nucleons are 
not just soft-mush but contain point-like constituents 
(now understood to be quarks) which are responsible 
for a power-law tail in high p± exclusive scattering, 
and in single particle inclusive production at high p±. 
As the physics of this phenomenon was explored ex- 
perimentally and theoretically, one of the most puzzling 
discoveries was the "Cronin effect" in which the cross- 
section for large-momentum transfer inclusive scattering 
was found to increase strongly with the atomic number 
of the target and the momentum transfer, compared to 
the behavior expected from the overall A 2 / 3 rescaling 
of the total cross section compared to p — p scattering. 
The dependence on A and p± can be parameterized as 



d 3 p j 



d 3 p 



(1) 



Fig. [T] shows a as a function of p±. The measurements 
were carried out at approximately 90° in the CM, and 
done for beam energies of 200, 300 and 400 GeV 
and targets ranging up to Tungston. Presumably, the 
Cronin effect results from multiple scattering, but it 
has never been satisfactorily explained theoretically 
and has not been further explored experimentally. The 
dependence of a on the quark-content of the secondary 
is particularly fascinating. 
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Fig. 2: The predicted normalized distributions of E±, for 
pions and kaons, of Pythia and QGSJet II, where E± is 
defined as yj m 2 + p\. The distributions are averaged 
over one million collisions of a lOOGeV proton upon a 
stationary proton. The top figure is the E± for secondary 
particles with a rapidity less than one. The middle 
figure is the E± distribution for secondary particles with 
rapidity between 1 and 2, and the bottom figure is the 
distribution for secondary particles with rapidity greater 
than 2. 



III. Difference between QGSJet II and Pythia 
FORp ± > 0.5 MeV 

The various event generators have different predic- 
tions for the production of particles with large p±. In 
particlular, Pythia, which is designed and tuned to de- 
scribe high p± particles, predicts a larger signal at high 
p± than QGSJet II. The difference in the predictions 
of the two event generators are illustrated in Fig. |2]for 
collisions of a lOOGeV primary proton on a stationary 
proton. We see that QGSJetll does not produce enough 
high p± particles. 

IV. Method 

Since current simulations do not include the "Cronin 
effect" and QGSJet does not produce enough high E± 
events even in p-p collisions, more realistic simulations 
would produce more secondaries at large scattering an- 
gle. A muon originating from a typical production height 
of 5 km, with a rs 10° production angle would arrive at 
1 km. To estimate the impact of increasing the number 
of high angle muons, we re-weight the secondaries in 
a suite of QGSJet simulations, by the factors shown 
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Fig. 3: The reweighting factors, as a function of p±, nec- 
essary to correct a target mass dependence of a = | to 
that of [1 1, averaged over the secondary types predicted 
by QGSJetll. 
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Fig. 4: The ratio of the E± distributions of Pythia and 
QGSJet II, using the same interactions as Fig. [2] for three 
bins in rapidity. The solid line covers the region in which 
we modified the secondaries in QGSJet II interactions 
to match the E± distribution on Pythia. 

in Figs. |3] and [4] denoting these weight increases by 

ronin 

(p±) and H Pythia (Ej_) respectively. 

Hcronin(p±) is found by taking the average of 
14. 4"% where i runs over particle species, weighted by 
the species abundance at high p±. Hp yt hi a (E±) is found 
by taking the ratios the the E± distributions for Pythia 
and QGSJetll. The E±_ distributions were generated for 
incident energies of lOOGeV, ITeV, and lOTeV. 

For every secondary particle produced in a QGSJetll 
interaction, the appropriate H is applied to the individual 
weight of the secondary particle. In order to conserve 
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Fig. 5: The ratio of the LDF with the Cronin- 
reweightings to the LDF when no change is made. The 
top figure is for a vertical shower, and the botton figure 
is for a shower inclined at 45°. For the muon ratios, 
what is plotted is the number density, while for photons 
what is plotted is the energy density. 

energy, the weight of all secondaries produced in the 
collision is decreased by a common factor chosen so 
the final energy equals the initial energy. This procedure 
effectively increases the number of high p± particles at 
the expense of multiplicity. 
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Fig. 6: The ratio of the LDF with the Pythia-reweightings 
compared to the LDF from QGSJet. The top figure is for 
a vertical shower, and the botton figure is for a shower 
inclined at 45°. For the muon ratios, what is plotted is 
the number density, while for photons what is plotted is 
the energy density. 



V. Results and Conclusions 



Fig. 7: The ratio of the LDF energy density of muons 
for both the Pythia-reweightings and the implementation 
of the Cronin effect, for a shower inclined at 45°. 

The impact on the Lateral Distribution Function, of 
the Cronin effect and of more accurately describing high 
pt particle production in the basic p-p collision, are 
shown in Figs. [5] and [6] Also shown in Fig.|7]are the ratio 
of the energy density of muons in the showers with and 
without reweighting. Evidently, the effect on the energy 
at large core radius of extra high-pt pions and the muons 
from their decay can be strong, but the number of high p t 
secondaries is probably simply too small to explain the 
observed signal in the surface detector. Thus, it seems 
that for the purposes of simulating extensive air showers, 
high-p t physics can safely be modeled crudely, without 
significant impact on the predicted surface array signal. 

Possibly more significant for the LDF is the produc- 
tion of relatively low energy pions by secondary particles 
from the primary collision which interact within the 
same nucleus before they escape. The results of a study 
of this effect will be presented at ICRC. 
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